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FOREWORD 


This  report  compares  the  performance  of  several  source  location 
algorithms:  Classical,  Adaptive,  Maximum  Entropy,  and  Eigenvector. 

The  results  should  be  useful  to  those  processing  sonar  signals  from 
uniform  line  arrays. 
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CHAPTER  1 
INTRODUCTION 


This  report  compares  the  resolution  and  the  detection  threshold 
of  various  multiple  source  location  algorithms.  Such  algorithms  are 
digital  signal  processing  techniques  that  use  signals  received  by  an 
array  of  sensors  and  provide  estimates  of  the  bearings  of  the  signal 
sources.  -We  will  confine  our  attention  to  linear  arrays  consisting 
of  omni-directional  acoustic  sensors.  Further  we  will  assume  the 
signals. to  be  narrowband. 

The  methods  under  consideration  are: 

1 )  Classical  method 

2)  Adaptive  method 

3)  Maximum  entropy  method  (approximate) 

4)  Eigenvector  method 

The  first  three  methods  involve  estimating  the  signal  power 
incident  on  the  array  from  various  bearings.  The  second  and  third 
method  involve  computation  of  the  maximum-likelihood  and  maximum- 
entropy  estimates  of  the  power  respectively  (reference  1  );  The 
eigenvector  method  is  based  on  certain  properties  of  the  eigenvectors 
of  the  correlation  matrix  of  an  array  illuminated  by  a  field  of 
discrete  sources  (references  2  and  3  ).  It  represents  the 
acoustic  field  as  a  finite  sum  of  plane  waves.— 

We  will  describe  these  methods  in  Chapter  2.  The  algorithms 
have  been  implemented  by  computer  programs  written  in  FORTRAN  and 
have  been  tested  with  simulated  data  on  the  VAX-11/780  system  at  the 
Naval  Surface  Weapons  Center,  White  Oak,  Maryland.  The  specific 
details  of  the  computer  programs  and  the  results  of  the  simulation 
are  described  In  Chapters  3  and  4  respectively.  In  Chapter  4  we  also 
compare  the  resolution,  detection  threshold,  and  computaional  effort 
of  the  four  methods. 
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CHAPTER  2 


DESCRIPTION  OF  THE  METHODS 


Before  proceeding  to  describe  the  multiple  source  location 
algorithms,  let  us  formulate  the  problem  under  consideration. 


Consider  a  linear  array  with  N  hydrophones.  The  array  responds 
to  signals  of  frequency,  f,  i.e.,  it  is  a  narrowband  array  with 
center  frequency,  f.  There  are  M  signal  sources  located  at  bearings 
angles  Oj  ,  02  ,..,  0M  with  respect  to  the  array  axis  (Figure  1)  and 
signal  powers  p  ,  p2,...,pM.  The  sources  are  assumed  to  be 
uncorrelated  and  the  incident  waves  are  planewaves.  The  noise  at  the 
hydrophones  is  assumed  to  be  isotropic  in  the  horizontal  plane  that 
contains  the  array  sensors. 


FIGURE  1.  GEOMETRY  OF  THE  PROBLEM 


ARRAY  AXIS 


Let  V(o  )  denote  the  steering  vector  corresponding  to  the  i-th 
source  of  planewave.  Then  V(0i)  is  given  by: 


V(0.) 


exp(j2  irdjCos  0/A) 
exp(.j2TTd2cosOi/A) 


exp(j2irdN_1cos01/A) 


(i) 


where  a  is  the  wavelength  of  the  incident  wave  and  0,di  ,d2,  •*»dN_i 
denote  the  hydrophone  positions.  We  denote  by  QN  the  NxN  covariance 
matrix  of  the  noise  signals  received  at  the  hydrophones.  The 
covariance  matrix  of  the  total  signal  received  at  the  array  elements 
Is  denoted  by  Q.  Since  the  sources  are  uncorrelated,  we  can  write  Q 
as 
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M 

Q  “  %  +  Z  Pi  v(V  V*  (G1)  (2) 

j  =  l 


where  *  denotes  complex  conjugate  tranpose. 

The  problem  of  Interest  to  us  can  now  be  stated  as  follows: 

given  Q,  compute  6  ,  1  =  1,2 . M,  the  bearing  of  the  sources.  We 

assume  that  M<N-1. 

Given  the  several  existing  methods  of  solving  this  problem,  we 
wish  to  compare  the  resolution  and  the  detection  threshold  of  these 
methods.  The  resolution  of  a  method  tells  us  how  close  two  sources 
can  be  before  they  become  Indistinguishable.  The  detection  threshold 
indicates  the  signal  power  level  above  which  the  source  can  be 
located. 

We  will  describe  four  different  methods  of  estimating  source 
locations  in  the  rest  of  this  chapter. 


CLASSICAL  METHOD 

The  estimate  of  power  incident  on  the  array  from  the  bearing  is 
taken  to  be 

P(0)  =  V*  (  0)QV  (  0)  (3) 

The  sources  can  be  identified  from  a  plot  of  P ( 0 )  against  6, 
Of  10°,  L80°|.  A  typical  plot  for  two  distinct  sources  Is  shown  in  Figure 
2.  The  sources  are  located  at  45°  and  120°. 


ADAPTIVE  METHOD 

In  this  method  we  use  the  max i mum- 1 i ke 1 i hood  estimate  of  the 
incident  power  distribution.  This  estimate  is  obtained  by  minimizing 
the  mean  square  error  for  unit  amplitude  planewave  signals  arriving 
at  the  bearing  for  which  the  incident  power  is  desired  (reference 
1).  The  adaptive  power  estimate  versus  bearing  is  given  by 
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i‘0>)  =  - — - -  (4) 

V*(6)  C!  V(0) 

A  plot  of  P ( 6 )  vs.O  should  reveal  the  source  locations.  In  our 
study, we  found  that  the  power  estimate  as  given  In  Equation  (4)  often 
fails.  Some  modifications  are  necessary  to  make  this  method  work. 
These  are  discussed  In  Chapter  3. 


MAXIMUM  ENTROPY  METHOD 


The  power  estimate  used  in  this  method  is  computed  using  the 
criterion  that  the  estimate  must  be  most  random  or  have  the  maximum 
entropy  of  any  estimate  which  is  consistent  with  the  measured 
covariance,  0  (reference  1  ).  The  power  distribution,  P(0),  is 
q  1  ven  by 

l>(0)  =  v*(0)  A  V(0)  (5) 

where  A  is  obtained  from  the  self-consistency  condition. 


1  r  V(0)V*(8) 
tt  J  V*(0)  A  V(0) 


d0  =  Q 


(6) 


In  general  no  method  is  known  for  solving  Equation  (6)  for  A.  A 
technique  due  to  Burg  (reference  T  )  can  be  used  to  solve  for  A  from 
Equation  (6)  under  the  following  conditions.  The  linear  array 
consists  of  equally  spaced  hydrophones  with  spacing  d=0.5x.  The 
covariance  matrix  Q  Is  Toeplitz  (we.  can  make  Q  Toeplitz  by  replacing 
each  element  In  a  sub-diagonal  by  the  average  of  all  elements  in  that 
sub -di a gona 1  ) .  In  this  case,  we  can  write  A  as 


A  =  YY* 


(7) 


Then  it  can  be  shown  that  reference  1 

i'  -  <rl  c  (8) 


where 

"  I 

0 

0 

e  -  0 


(9) 


LoJ 

We  should  note  that  A  can  be  expressed  in  the  form  given  by  Equation 
(7)  whenever  the  line  array  is  composed  of  equally  spaced  sensors. 
For  an  equispaced  line  array  with  d  >  0.5A  we  get  ambiguous  estimates 
for  the  source  bearing. 
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There  are  no  known  results  for  the  case  d  <  0.5A.  We  can, 
however,  use  an  approximate  method  in  this  case.  This  method  is 
described  in  Chapter  3. 


EIGENVECTOR  METHOD 

Here  we  determine  the  source  locations  from  certain  properties 
of  the  eigenvectors  of  the  covariance  matrix,  Q,  (references  2  and 
3  ).  To  start  with,  suppose  the  noise  signals  at  the  array  elements 
are  uncorrelated  with  identical  variance,  a2.  Then  Q  =  a2I,  where  I 
is  an  NxN  identity  matrix.  Then; 

M 

qn  =  02I  +  £  Pi  v(e.)v*(0.)  (10) 

The  algorithm  is  based  on  the  following  theorem  due  to  C  '•"t  oni  and 
Godara  (reference  3  ). 

Theorem:  Let  W  be  the  eigenvector  of  Q  corresponding  tc  -:e  smallest 

eigenvalue.  Then  if  M^N-i 

W*V(oi  )  =  0,  1  =  1,  2  , . . .  ,M  (11) 

We  note  that  this  theorem  can  be  applied  if  the  sources  are 
correlated.  In  any  case,  the  theorem  suggests  an  algoritnm  for 
finding  the  source  locations  from  the  covariance  matrix  Q.  The 
algorithm  is  as  follows: 

1)  Find  the  eigenvalues  and  eigenvectors  of  Q 

2)  Find  W,  the  eigenvector  of  Q  corresponding  to  the 
smallest  eigenvalue. 

3)  Solve  for  the  roots  of 
1 

exp(  jZirdjCose/x) 
exp(  j2ird2coso/x) 

•  • 

exp(j2  dN _-|Cose/A) 

These  roots  give  the  source  bearings, 

4)  Put  the  V(e1)  into  Equation  (10)  and  find  power, 

Plt  of  each  source.  Then  those  sources 
corresponding  to  power  greater  than  some 
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1 


threshold  value  can  be  considered  to  he  actually 
present. 

The  above  algorithm  Is  based  on  the  assumption  that  the  noise 
signals  are  uncorrelated  (I.e.,  Q  Is  proportional  to  I).  We  can  use 
the  same  algorithm  for  correlated  noise  signals  provided  we  p re whiten 
the  covariance  matrix,  Q,  and  modify  the  eigenvector,  W,  before 
solving  Equation  (12).  These  modi f i ca t i ons  are  indicated  below: 

1)  Since  Q  is  Hermitian,  we  can  write  it  as 
Q  =LL* 

2)  Modify  Q  by  transforming  it  as  follows: 

Q  L  Q  ( L  )* 

3)  Find  W,  the  eigenvector  corresponding  to  the 
eigenvalue  of  the  transformed  Q.  Now  transform  W  as 
follows: 

W*  <-  W*L 

and  use  this  W*  in  Equation  (12). 

The  implemented  version  of  this  algorithm  turned  out  to  be 
numerically  unstable  (Chapter  3).  A  modified  algorithm  (reference 
4  )  was  used  in  the  final  implementation.  We  will  describe  this  in 
the  following  chapter. 


CHAPTER  3 


IMPLEMENTATION  OF  THE  ALGORITHMS 


In  Chapter  2  we  have  described  four  methods  of  estimating  source 
hearings  using  a  linear  hydrophone  array.  In  our  simulation  we 
allowed  the  sensor  positions  to  be  perturbed  by  a  small  Gaussian 
random  displacement  from  their  nominal  positions.  We  use  spacings 
of  ,3A  (\  is  the  wave  length  of  the  incident  wave).  This  in  turn 
gives  us  a  rather  singular  covariance  matrix.  Modifications  are 
therefore  needed  for  implementation  of  the  algorithms. 


MODIFICATIONS  TO  THE  CLASSICAL  METHOD 

Classical  method  has  poor  resolution  and  detection  threshold. 
To  improve  the  detection  threshold  we  do  the  following: 

1)  estimate  s  (power  of  the  background  noise) 

2 )  compute  R=  Q  -  sQ  , 

N 

3)  then  P( o)  =  V*( 0)RV(e ) 


MODIFICATIONS  TO  THE  ADAPTIVE  METHOD 

When  implemented  using  Equation  (4),  the  power  distribution 
versus  bearing,  P(o),  given  by  the  adaptive  method  is  erroneous. 
This  is  due  to  numerical  instability  in  inverting  the  covariance 
matrix,  Q.  Q  is  ill-conditioned  in  case  of  small  spacings  (over 
sampled),  i.e.,  the  spread  between  its  largest  eigenvalue  and  its 
smallest  one  is  very  large.  The  inevitable  uncertainties  in  the 
hydrophone  locations  produce  large  percentage  changes  in  the  smaller 
eigenvalues  of  Q.  This  in  turn  affects  the  evaluation  of  P(e).  Our 
experience  has  shown  that  the  degrading  effects  of  this  problem  can 
be  alleviated  if  we  modify  the  diagonal  elements  of  Q  as  follows: 
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MODIFICATIONS  TO  THE  MAXIMUM  ENTROPY  METHOD 


For  an  equally  spaced  array  with  spacing  d<  0.5A  and  for  a 
Toeplltz  covariance  matrix  we  use  the  following  method  to  find  the 
source  bearings. 


Let  a  be  defined  by  Equation  (7).  Then  it  can  be 
(see  Appendix  A) 


shown  that 


r  ft  g  ‘n 


(13) 


where 


r-'«(,o 


i 

(> 


(,<> 


(14) 


L  I  ((n-i  )i-;  )J 


2nd 

and  J  is  the  zeroth  order  Bessel  function  of  the  first  kind. 

t  I 

The  maximum  entropy  estimate  of  the  source  bearings  is 
given  by  the  angles  where  the  local  minima  of  | v* ( 3 ) r  |  occur. 


(15) 


then 


MODIFICATIONS  TO  THE  EIGENVECTOR  METHOD 

The  eigenvector  method  in  the  modified  form  is  as  follows:  Let 

n  denote  the  covariance  matrix  of  an  assumed  noise  model,  then 

N  M 

Q  =  sQn  +  £  P1V(0.)V*(O.)  (16) 

i=l 

where  s  Indicates  the  noise  level.  Let  Q  be  Toeplitz.  The  following 
algorithm  gives  us  an  estimate  of  the  source  bearings. 

1)  Find  s,  an  estimate  of  s  and  write 

R  “  0  -  ?Qn  (17) 

2)  Find  eigenvectors  and  eigenvalues  of  R 

3)  Estimate  M,  the  number  of  sources  present.  This 
can  be  done  by  studying  the  eigenvalues  of  R. 

4)  Find  the  source  bearings  by  solving: 

min  |  Y  *V ( 0 )  |  (18) 

n 

where  y  is  an  Nx(N-M)  matrix  formed  by  (N-M)  eigenvectors  of  R 
corresponding  to  its  (N-M)  smallest  eigenvalues  (reference  4  ). 
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1)1  f  A  1 1  S  OK  I  MP1  K  Ml  NT  A  F  I  ON 

The  soutch  location  algorithms  are  tested  using  simulated  data. 
The  usual  approach  is  to  generate  sequences  of  data  vectors  for  the 
linear  array  and  then  compute  the  covariance  matrix  Q  from: 

K 

•>  -  l  £  x(i)x*(i)  (19) 


where  X(i)  denotes  the  i-th  data  vector.  In  the  simulation,  we  ncve 
opted  to  generate  Q  directly.  Our  method  incorporates  the  effects  of 
computing  Q  using  Equation  (19)  with  finite  K.  (See  Appendix  B.) 

The  noise  covariance  matrix,  QN,  is  generated  by 


o 

'N 


l 

N 


k 

y  v(ojv*(u.) 

i-i  1 


(20) 


where  V(  u.)  is  given  by  Equation  (1)  with  u-  =  ^  (i-0.5).  A 
deterministic  0  is  then  computed  from  Equation  (2),  with  P.  and  6i  , 
i  1, M  appropriate  for  the  simulation  being  made.  In  order  to 
introduce  the  statistical  effects  of  computing  Q  using  Equation  (19) 
with  finite  K,  we  use  the  following  algorithm. 


1)  factor  Q  as  :  Q  =  UU*  where  U  is  the  Choleski  factor 
of  0 

2)  generate  a  perturbed  identity  matrix,  I,  (see 
Appendi x  B ) 

3)  Define  Q  as  :  Q  =  UI'O* 


4)  To  smooth  Q  we  average  all  duplicated  inter-element 
correlations  in  the  covariance  matrix.  In  the  special 
case  of  equal  spacing  this  involves  averaging  the  elements 
down  each  diagonal  resulting  in  a  Toeplitz  matrix.  This  is 
necessary  to  obtain  a  covariance  estimate  consistent  with 
the  model  we  use. 
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CHAPTER  4 

RESULTS  OF  THE  SIMULATION 


In  order  to  compare  the  performance  of  the  four  algorithms,  we 
generate  a  sequence  of  covariance  matrices  corresponding  to  a  number 
of  time  cuts.  These  covariance  matrices  contain  a  target  track 
toqether  with  background  noise.  In  generating  these  covariance 
matrices  we  assume  that  the  position  of  each  array  sensor  has  been 
perturbed  by  a  small  Gaussian  random  displacement  (o  =.01  on  both 
coordinates).  In  our  simulation  we  use  a  linear  array  with  8 
sensors  of  spacings  equal  to  ,3y.  Each  of  the  implemented  algorithms 
is  applied  to  this  data  and  the  estimated  target  track  is  plotted. 
The  detection  threshold  and  the  re solution  of  an  algorithm  is  found 
from  the  target  track  computed  by  that  algorithm.  We  vary  the  signal 
power  level  to  study  the  resolution.  We  also  vary  the  number  of  post 
i n tegrat i ons ( 1 . e . ,  the  number  of  averaging)  and  study  their  effect  on 
the  resolution  and  detection  threshold  of  the  algorithm. 


DETECTION  THRESHOLD 

The  detection  threshold  is  defined  as  the  signal  power  level 
above  which  the  source  can  be  located.  In  order  to  compare  the 
detection  threshold  of  the  algorithms  it  is  necessary  to  calculate, 
for  each  method,  another  threshold  which  we  will  call  the  "decision 
threshold".  An  algorithm  would  indicate  the  presence  of  a  source  at 
a  certain  bearing  if  the  power  from  that  direction  exceeds  the 
"decision  threshold"  of  that  algorithm. 

The  decision  threshold  of  each  algorithm  is  found  in  the 
following  manner.  We  generate  a  large  number  (**1  ,000)  of  covariance 
matrices  containing  noise  only.  For  each  covariance  matrix,  an 
algorithm  would  compute  the  peak  power.  The  maximum  in  the  sequence 
of  peak  power  is  used  as  the  decision  threshold  of  the  algorithm  for 
a  false  alarm  rate  of  approximately  .001. 

To  find  the  detection  threshold  of  an  algorithm,  we  generate  a 
sequence  of  covariance  matrices  containing  a  sawtooth  track  (see 
Figure  3).  At  intervals  of  20  time  cuts  the  signal  power  level  is 
decreased  by  1  dB.  An  algorithm  is  then  applied  to  the  sequence  of 
covariance  matrices.  At  each  time  cut,  the  algorithm  produces  a 
source  location  provided  the  power  received  from  that  direction 
exceeds  the  decision  threshold  of  that  algorithm.  The  detection 
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threshold  is  given  by  the  signal  power  level  at  which  the  algoritnir. 
finds  the  source  approximately  fifty  percent  of  the  time. 

TABIC  1.  COMPARISON  OF  CLASSICAL,  ADAPTIVE,  MAXIMUM  ENTROPY, 

EIGENVECTOR  FOR  DETECTION  THRESHOLD 


ll  of  post 
Intep.rot  ionS 

Classical 

Adaptive 

M.E. 

Eigenvector  ' 

1 

1,000 

-13.5  dB 

More  than 

-11  dB 

-li  dB 

-H-5  qB 

100 

-  8.5  dB 

More  than 
-  7  dB 

-  7 dB 

-  7.5  dB 

i 

Table  1  indicates  that  the  classical  method  performs  better  than 
the  others  in  detecting  a  single  signal.  The  adaptive  method  has 
poor  detection  threshol d, and  it  gives  a  biased  result.  We  also 
learned  that  the  perturbation  of  sensor  positions  does  not  change 
performance  significantly.  However,  the  variance  of  estimated 
bearings  increases  due  to  perturbation  of  array  sensor  positions. 


RESOLUTION 

The  resolution  of  each  algorithm  is  found  as  follows:  a 

sequence  of  covariance  matrices  is  generated  which  contain  an 
Interference  line  and  a  track  approaching  this  line.  The  separation 
at  which  the  two  sources  become  barely  distinguishable  gives  the 
resolution  of  the  algorithm.  At  each  time  cut,  a  reasonable 
quantitative  measure  of  the  resolvability  would  be  minimum  angle 
separation  above  which  two  signals  are  resolved  and  below  which  they 
are  not  (i.e.,  minimum  angle  between  two  signals  such  that  they 
resolve  in  the  sense  that  the  estimated  power  spectra  display  two 
distinct  peaks  for  classical,  adaptive  and  maximum  entropy 
algorithm). 

The  resolution  of  the  classical  method  can  be  estimated  as 
follows:  Suppose  we  have  two  signals,  one  located  at  x0  and  the  other 
at  -x  .  Let  V(x)  be  the  steering  vector  at  x  where  x=cos(fi)  (i.e., 
the  x-th  element  of  V ( x ) =  exp  ( i  2tts  xk/X  )  where  s  is  the  separation 
between  array  sensors  and  X  denotes  the  wavelength  of  the 
incident  wave). 
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Then  the  power  spectrum  at  x,  n(x),  is  given  by 


. ■  ( * )  v* ( x )  |  v ( x  ) v*(x  )  +  v (-  x  )  |  V(x) 

o  o  o 


from  which  we  can  write 


o(o)  =  2 


J"  sin(27TsxoN/A) 


sin(7rsx 


oW/»' 

o«>  . 


where  N  is  the  number  of  array  sensors  and 

[sin(2irsxoN/A)  "1 

7in(27Tsx0o/A)  J 


Let  6  =  ttsx  N/A 
o 

then  - jr  o(o)«s  2(sin  0/8)2  (since  sin(nsx  /A)^3  for  N»l) 
N  o 

and  ~N  (xo)*l  +  [sin28/28] 2 


therefore 


o(o) 


402  +  (sin20) 2  1 

8sin20  2 


K 


sin0 


)2 


+  cos 2  8 ] 


We  have  resolution  when  o(x  )/o(l)>l  or  minimum  resolution  occurs  when 
o(x„)/o(o)  =  1,  in  this  case0,  we  get:  8  =  .9^  .  Minimum  resolution  for 
classical  sss  21°  for  N  =  8,  s/A  =  .3. 

To  study  the  effect  of  the  signal  power  on  the  resolution,  we 
have  Investigated  a  number  of  simulations  using  different  signal 
power.  The  main  result  which  we  have  expected  and  observed  is  that 
the  stronger  the  signal  power  the  better  the  resolution  (with  the 
exception  of  the  classical  algorithm). 
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r All L E  2.  COMPARISON  OF  CLASSICAL,  ADAPTIVE,  MAXIMUM  ENTROPY, 

EIGENVECTOR  FOR  RESOLUTION  USING  1000  POST  INTEGRATIONS 
( 5;  e  e  Figure  4  ) 


RESOLUTION  OF 


Signal  Power 

Classical 

Adaptive 

M.E. 

Eigenvector 

0  ill> 

21° 

15° 

9° 

4° 

-0  ill) 

21° 

19° 

12° 

8° 

-12  .11) 

21° 

23° 

21° 

11°  j 

TABLE  3.  COMPARISON  OF  CLASSICAL,  ADAPTIVE,  MAXIMUM  ENTROPY, 

EIGENVECTOR  FOR  RESOLUTION  USING  100  POST  INTEGRATIONS 

RESOLUTION  OF 


Signal  Power 

Classical 

Adaptive 

M.E. 

Eigenvector 

0  db 

22° 

15° 

0 

O 
«— 1 

7.5° 

-b  <ll) 

22° 

19° 

12° 

9.5° 

-9  dl> 

22° 

20° 

15° 

12.0° 

Tables  2  and  3  indicate  that  the  eigenvector  method  performs  consis¬ 
tently  better  than  the  others  in  resolving  targets.  Classical 
method  has  poor  resolution,  but  the  resolution  does  not  depend  on 
power  in  this  case.  Therefore  in  case  of  very  weak  signal  the 
classical  algorithm  can  be  used  successfully  while  the  others  are 
no  longer  applicable. 
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COMPUTATION  T  I  Mb 

We  have  implemented  the  four  source  location  algorithms,  viz, 
classical,  adaptive,  maximum  entropy  and  eigenvector  algorithm. 

The  programs  were  written  in  FORTRAN  with  no  particular  effort  made 
to  make  them  time  and  memory  efficient.  We  will  estimate  (Table  4) 
the  number  of  operations,  in  each  algorithm,  in  flops  (a  "flop"  is 
roughly  equivalent  to  the  amount  of  time  required  to  perform  one 
multiplication,  one  addition,  and  perhaps  a  few  address  calculations). 

TABLE  4.  COMPARISON  OF  CLASSICAL,  ADAPTIVE,  MAXIMUM  ENTROPY, 


EIGENVECTOR  FOR 

COMPUTATIONAL  TIME 

(FLOPS) 

Classical 

Adapti ve 

M.  E. 

Ei gen  vector 

15N 

1GN+3N2 

15N+3N2 

15N+N3 

N:  number  of  array  sensors. 


In  all  fairness,  we  should  make  it  clear  that  these  estimates  of 
number  of  computations  necessary  to  implement  each  algorithm  are  by 
no  means  exact  yet  the  numbers  do  indicate  that  the  eigenvector 
algorithm  takes  much  more  time  to  execute  than  the  other  algorithms. 
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CHAPTER  5 
CONCLUSIONS 


The  results  of  a  study  of  source  location  algorithms  have  been 
presented.  The  algorithms  were  compared  on  the  basis  of  their 
ability  to  resolve  closely  positioned  sources  and  locate  low  power 
siqnals.  The  computational  effort  involved  in  implementing  the 
algorithms  was  also  considered. 

The  classical  method  has  the  best  detection  threshold, 

approxi mately  1  to  3  db  less  than  that  of  the  others.  The 

eigenvector  method  comes  next,  followed  by  the  maximum  entropy  method 
and  lastly,  the  adaptive  method.  It  was  found  that  increasing  the 

number  of  post  integrations  from  100  to  1,000,  in  spite  of  coordinate 

perturbations,  reduces  the  detection  threshold  by  nearly  the  expected 
amounts. 

The  eigenvector  method  provides  the  best  resolution;  the  maximum 
entropy  and  the  adaptive  method  follow  with  gradually  worsening 
resolution.  The  classical  method  can  resolve  signals  only  when  they 
are  very  far  apart  (i.e.,  21°  apart  in  the  case  of  8  sensors  with 
0.3*  spaci ngs  ) . 

The  classical  method  requires  the  least  amount  of  computation 
and  is  very  cost  effective  provided  resolution  is  not  important.  The 
comput.  at  i  ona  1  efforts  involved  in  the  adaptive  and  the  maximum  entropy 
methods  are  about  the  same.  The  eigenvector  method  is  rather 
expensive  compared  to  the  other  methods. 

In  our  simulation  an  array  of  eight  sensors  was  used.  For  arrays 
involving  a  large  number  of  sensors,  our  conjecture  is  that  the  maximum 
entropy  methods  provides  the  best  performance  at  reasonable  computa¬ 
tion  cost.  Further  investigation  is  necessary  to  confirm  this 
conjecture. 
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APPENDIX  A 

DERIVATION  OF  APPROXIMATE  ME  METHOD 


In  this  appendix  we  derive  an  approximation  for  the  M.E.  method, 
provide  justifications  for  this  approximation  in  simple  cases. 

We  have: 

7T 

Q  =  f  do  V(0)  V*  (0)  o(0) 

-TT 

where  a ( O )  is  the  power  density  at  bearing  0 

Vk(0)  =  ei^cos0k  (k**1  element  of  V ( e ) ) ,  where  8  = 

then 

qk1  ■  /ei8(k-' )cosoa(o)  de 

-7f 

Let  x  =  0coso*  we  get: 


Let  f(x)  = 


}  -lx(>'li 

-8  V^2  ~  *2 

-.■-L=  x2  <  e2 
V8?-  x2 


=  o  x2  >  e2 


Express  f(x)  using  Fourier  series,  we  obtain 


inx 


f(x)  =  L  cne 

n=-<»> 

>  dx  e-in* 

cn  /  2, 

-71 

cn  «  J0(8n) 


=  }  do  e-i0ncose 


Then 


oki  -  z  y«o  /  e,x(k''  * n)  fix) 

n»-">  -n 


dx 


A-l 


We  also 


(A.  1 ) 


(A. 2) 
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From  equations  (6)  and  (7),  we  have 
1 


V*(»)y>*V(0) 

ii  i(k-l+n)x, 

Qkl  =  E  J0(Bn)  f  - — 

kl  n~-«  0  -Jn  V*(x)yy  *V(x) 

r  Wv 

Let  dm  =  J  (V*y)(y*V) 


Then  a  =  a  * 
-m  m 


Qkl  =  J  VB">  al-k-n 


n=-<» 


t  Sv,-k.n 


Let  0 


i=  .?  v, 


1=0 


i+1 


It  can  be. shown  that 
o  l<o 


°1 


1 


l=o 


For  1 >  o 


(Y  *)2 

y  O  2  ~  ^ 

(Y0*)2 


(y.*)2  -  Yn*  Y 


*  v  * 
0  J2 


Therefore, 


(y  *y 

'  'o 


N-  1 

E  Qk-|  Y  ^ 
1=0  Kl  1 


■rj.  [j0(Bk)  -(£)  *  J0(B(k 
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Taking  only  the  1st  term  of  the  series  and  ignoring  the  constant  O/y  Q*)  we 

get: 

“.(■)  1 

,  Jo"»  I 


|j0((N-l)6) 
Let  J  be  defined  as: 

"V°> 

J0(B) 

0  .  V*’ 


Lj0((n-i)b1J 

If  the  background  noise  is  isotropic  in  the  plane  which  contains  the 
array  sensors  then  J  is  the  first  column  of  the  covariance  of  background 
noise.  In  other  words, 

J  ■ 

where  is  the  covariance  of  background  noise  and 
—  — 

1 

0 

0 


L°J 

Suppose  we  have  no  signal  (i.e.,  Q  =  Qn)  then  with  the  approximation  in 
Lquation  (6) ,  we  get: 

v  =  Q_1J  =  QN]  QN  e  =  e 

Therefore  V*(z)  y  =  z^  ^ 

So  with  this  approximation  of  M.E.  we  obtain  a  flat  power  spectral  (i.e., 
all  poles  at  the  origin)  as  when  there  is  no  signal. 


Suppose  we  now  have  one  planewave,  then: 

Q  =  Qn  +  pV(0o)V*(0o) 

«''■  -  WCp  Qn1 
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where  A  =  V*( oq)  Q^V(oo) 

y  =  Q"V 

therefore: 

k  -  r+lp  <  v<Vv*<eo>  «N  e 

V  -e  -  Q-’v(e0)V*(60)  e 

P 

■ «  -  <T-f  >  *(60) 

p 

v*(0)  a  =  1  -  (y^-)V*(e)  q;1V(60) 

p 

The  second  term  in  the  right  hand  side  of  the  above  equation  is  the  optimum 
processing  for  one  plane  wave.  So  using  this  approximation  we  process  the  signal 
optimumly  in  case  of  one  plane  wave. 
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APPENDIX  U 

COVARIANCE  MATRIX  SIMUI.ATION 


In  this  Appendix  we  provide  justifications  for  generating  covariance  matrix 
In  our  a  I  mu  I  at  ion . 


J,oL  v(k)  be  the  k-tli  observation  vector 


Then : 

1  M 

Q  =  £  £  y00y*00 

k=l 

y(k)y*(k)  can  be  factored  as:  yy*  =  UU* 
bet  z(k)  be  such  that  z(k)z(k)*  =  I 
Then  y(k)  can  be  written  as: 

Y(k)  =  Uz (k) 


M 

£y(k)y*(k)  = 
k~  1 


r M 

1  52  z(k)z*(k) 

Lk=i  J 


u* 


=  M  52  y(R)v*(R)  =  u  ^  L 

L  ^=1 


M 

52 

k=l 


z (k) z*(k)  U* 


()  *  UiU* 

M 

where  T  =  -  52  z(k)z*00 

k=i 


(B .  1) 


(B.2) 


So  Equations  (1)  and  (2)  account  for  Steps  (1),  (2),  and  (3)  (in  Chapter  3)  in 
generating  covariance  matrices  with  statistics. 


To  generate  a  perturbed  identify  matrix,  l,  let  us  define  such  that 


I  1  * 

”M  r! 


■t 


z(k)z(k)* 


k=  I 


In  other  words,  1^  is  a  Choleski  factor  of  £  z(k)z(k)* 


k=l 


B-L 


NSWC  TR  83-12 


It  f.m  lx-  shown  tli.it  1^  is  a  lower  Triangular  matrix  where  diagonal  elements 
are  eh  I  d i st r I  billed ,  and  off  diagonal  elements  are  Rayleigh  distributed  with  unit 
variance.  All  the  elements  are  independent. 
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